Localized Magnetic States in Graphene 
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We examine the conditions necessary for the presence of localized magnetic moments on adatoms with inner 
shell electrons in graphene. We show that the low density of states at the Dirac point, and the anomalous 
broadening of the adatom electronic level, lead to the formation of magnetic moments for arbitrarily small local 
charging energy. As a result, we obtain an anomalous scaling of the boundary separating magnetic and non- 
magnetic states. We show that, unlike any other material, the formation of magnetic moments can be controlled 
by an electric field effect. 
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Graphene, a two-dimensional (2D) allotrope of carbon, has 
singular spectroscopic and transport properties d 0. H, 0] 
due to its unusual electronic excitations described in terms of 
massless, chiral, "relativistic" Dirac fermions [5]. Besides be- 
ing a possible test bed for relativistic quantum field theory [6], 
graphene has a great technological potential due to its struc- 
tural robustness, allowing extreme miniaturization [7], and a 
flexible electronic structure that can be controlled by an ap- 
plied perpendicular electric field JUl. 

In this paper we show that graphene has also potentiality for 
spintronics, that is, independent control of the charge and the 
spin of the charge carriers (19]. Unlike diluted magnetically 
semiconductors (DMS) ifioll where the location of the mag- 
netic ions is random and hence unpredictable, adatoms can 
be positioned in graphene using a scanning tunneling micro- 
scope (STM) 111 111 . Furthermore, as we are going to show, the 
magnetic properties of adatoms such as size of the magnetic 
moment and Curie temperatures can be controlled by an ex- 
ternal electric field, an effect unparalleled in condensed matter 
systems. 

The basic model for the study of magnetic moment forma- 
tion in metals is the well-known Anderson impurity model 
1 12J • In this model an ion with inner shell electrons with en- 
ergy e hybridizes, via a hopping term of energy V, with a 
conduction sea of electrons. While the conduction electrons 
are described by a Fermi liquid with featureless, essentially 
constant, density of states (DOS), the impurity ion is assumed 
to be strongly interacting. The Coulomb energy required for 
double occupancy of an energy level in the ion is given by U. 
Anderson showed that when eo is below the Fermi energy, //, 
and the energy of the doubly occupied states, eo + U, is larger 
than /i, a magnetic state is possible if U is sufficiently large 
and/or V sufficiently small. 

Here we apply the Anderson model to graphene and show 
that the energy dependence of the DOS leads to anomalous 
broadening of the adatom level and strongly favors the forma- 
tion of local magnetic moments. In particular we show that, 
unlike the case of ordinary metals, this anomalous broaden- 
ing allows the formation of magnetic states even when eo is 
above the Fermi energy at relatively small U . We also find, in 
contrast with the usual metallic case, an anomalous scaling of 



the magnetic boundary separating magnetic and non-magnetic 
impurity states. Finally, we establish that the local magnetic 
moments can be mastered by the application of an external 
gate voltage, leading to a complete control of the magnetic 
properties of adatoms in graphene. 

We consider an impurity atom adsorbed on the surface of 
the graphene sheet, on top of a carbon (see Fig. [T]). The tight- 
binding Hamiltonian of the electrons in graphene is 

H TB = ~tJ2Yl [4(Ri)k(Btf) +H.c] , (1) 

where <z CT (R,) (6 CT (R,)) annihilates and electron with spin 
a =1, i on sublattice A (B) at position R^, stands for 
summation over nearest neighbors, and t (« 2.7 eV) is the 
nearest neighbor hopping energy. In momentum space, we 
have (we use units such that h = 1): 
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where <j>(k) = £ f e ik? , with <fi = a(x/2 + \/3/2y), 5 2 = 
a(x/2 — y/3/2y) and 5^ = —ax are the nearest neigh- 




Figure 1: (color on line) (a) Honeycomb lattice with an impurity 
atom. Black: sublattice A; White: sublattice B. Intersection of the 
Dirac cone spectrum, e±(k) = ±Vpk, with the localized level spec- 
trum, E f (k) = eo: b) e > 0; (c) e < 0. 
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bor vectors. Diagonalization of the Hamiltonian (fJJ gener- 
ates two bands, e±(k) = ±t|</>(k)|, which can be linearized 
around the Dirac points K at the corners of the Brillouin zone: 
e±(K + q) » ±i>p|q|, where Vp = 3ta/2 (ss 10 6 m/s) is the 
Fermi velocity of the Dirac electrons. 

The hybridization with the localized orbital of the impurity 
atom on a given site, say, on sublattice B, is given by: Hy = 

VY,Ma b *{°) + H - C A > where f° (/<r) annihilates (creates) 
and electron with spin a =f, J, at the impurity. In momentum 
space we have: 



(3) 



where iVj, is the number of sites on sublattice B contained in 
the expanded unit cell of graphene with the impurity. 

The Hamiltonian of the localized orbital is described by a 
single level, Hf = eo J2 a ftf<?- The electronic correlations in 
the inner shell states can be described by a Hubbard-like term: 
Hjj — Uftf^flfi. Following Anderson, we use a mean- 
field decoupling of the interaction, Hjj — > Un- a f^f a — 
Un^ni, where n a = (flf a ) is the occupation for each of the 
two spin states. The Hubbard term can be absorbed into the 
definition of the local impurity energy, Hf = e <jflf<r i 
where e a = eo + Un- a is the energy of the localized electrons 
in a given spin state in the presence of a charging energy U. 

The formation of a magnetic moment is determined by the 
occupation of the two spin states at the impurity, n a . A lo- 
calized moment forms whenever n-f 7^ Uj The determination 
of n a requires the self-consistent calculation of the density of 
states at the impurity level, Pff(u>), which incorporates the 
broadening of the impurity level due to hybridization with the 
bath of electrons in graphene. The occupation of the impurity 
level is given by: 



(4) 



The Green's function of /-electrons is: Gff. a (t) - 
-i(T [U(t)fl(0)]), and its retarded part can be written as: 



Gf/,»= [ W -e <r -E«( W )+iO+] \ (5) 



where 



E*M = (yV^)£G°*(p, 



(6) 



is the self-energy of the /-electrons, which is defined 
in terms of the non-interacting Green's function of the 
graphene electrons in a given sublattice, G® bcr (p,t) = 
-z(T[& ffp (t)&t p (0)]> : 

G&*(p, u>) = lo/(lo 2 - v 2 F \p\ 2 + zO+signM) . (7) 
In this case (O becomes: 



(8) 



where D is a high-energy cut-off of the order of the graphene 
bandwidth (D w 7 eV). We choose the cut-off D using the 
Debye prescription, i.e., conservation of the number of states 
in the Brillouin zone after linearization of the spectrum around 
the K point. We assume <C D, where band effects related 
to the exact definition of the cut-off are not important. 

The real part of defines the quasiparticle residue 

Z-^uj) = 1 + {V 2 /D 2 )\n(\D 2 -cj 2 \/uj 2 ) of the /- 
electrons, while the imaginary part gives the broadening of 
the localized level due to the hybridization. As expected, the 
anomalous character of the problem is explicitly manifested 
in the linear dependence of the broadening with the energy, 
which is proportional to the electronic DOS in graphene. Fur- 
thermore, notice that Z(oj) vanishes at ui — > 0. Replacing Eq. 
© into Eq. @ gives the density of states of the localized 
level, pff, a (v) = ~ 1 / 7rImG '//. C r( w ) : 



1 AM0(£>-M) 



Tr[Z^(uj)uj-e a ] 2 + A 2 lo 2 



(9) 



where A = ttV 2 /D 2 is the dimensionless hybridization. 

Notice that, unlike the case of impurities in metals, the im- 
purity density of states is not a simple lorentzian. The im- 
purity DOS, @, is peaked around the quasiparticle pole at 
€ a = and u> — ► 0. We can expand Z{u) around the sin- 
gularity at uiqZ^^-^q) ~ e CT for uiq — > 0, where we may 
approximate Z(loq) ~ Z(e a ) except for a double-logarithmic 
corrections that can be safely ignored lfTil . The anomalous 
broadening gives rise to a logarithmic divergence in the ultra- 
violet when the DOS of the level is integrated in ©, 
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where Z(e a ) = Z a and p a = pZ a 1 . The term G CT contains 
the contribution coming from the cut-off regularization: 



where 7 = sign(/i), and 



Q a = Z a — In 
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E a = V /(e ff -^ 1 )2 + M 2 A 2 
W a = J{DZ^+e a 
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In Fig. 2 we show the boundary between magnetic and 
non-magnetic impurity states as a function of the parameters 
x = DA/U and y — (p — eo)/U. Notice that the magnetic 
boundary is not symmetric between the cases where the im- 
purity is above (eo > 0) or below (eo < 0) the Dirac point. 
Moreover, unlike the metallic problemJT^] the boundary is not 
symmetric around y = 0.5. This reflects the particle-hole 
symmetry breaking due to the presence of the localized level. 
In the case where eo > (see Fig. 2a), the magnetic boundary 
crosses the line y — 0, and the level magnetizes even when 
the impurity is above the Fermi energy. This is understood by 
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Figure 2: Boundary between magnetic and non-magnetic impurity 
states in the scaling variables x and y for eo > (a) and eo < (b). 
Circles: |e |/-D = 0.029, V/D = 0.14; Squares: e /D = 0.043 
and V/D = 0.14; Triangles:|e |/D = 0.029, V/D = 0.03. The 
upturn close to y = 1 and x — > on panel b) is not visible in this 
scale when V is very small (triangles). See details in the text. 



the fact that the hybridization leads to a large broadening of 
the impurity level density of states (with a tail that decays like 
1 /lo) that crosses the Fermi energy even when the bare level 
energy is above it. In the opposite case of eo < a similar ef- 
fect occurs with the crossing of the magnetic boundary along 
the y = 1 line, something that also does not occur in ordi- 
nary metals 114f| . This implies even when the energy of the 



doubly occupied state is below the Fermi level (eo + U < fi), 
because of the large broadening, the impurity magnetizes if 
U is not too large or too small. The up turn close to y — 1 
and x <C 1 in the eo < case only reflects that in this limit ( 
U, [i 3> eo, for finite eo) the physics of the Dirac points is ir- 
relevant and we recover the usual Anderson model in ordinary 
metals, where the transition curve approaches the point x = 0, 
y = 1 from below lll2ll . This picture is physically consistent 
with a renormalization group calculation at fixed /i = 115 ] . 

The dependence of the scaling of the magnetic boundary 
with eo and A (see Fig. 2) shows that the size of the mag- 
netic region grows as |eo| approaches the energy of the Dirac 
points. In this situation the DOS around the localized level 
is suppressed, favoring the formation of a local magnetic mo- 
ment. In particular, in the limit where the level is nearly at the 
Dirac point (|eo| — > 0), the level nearly decouples from the 
bath and the impurity can magnetize in principle for any small 
finite charging energy U. On the other hand, the magnetic re- 
gion shrinks in the y direction as the hybridization parameter 
A grows (see Fig. 2). In the limit of A — * and U finite a 
local magnetic moment forms whenever < y < 1, as in the 
case of an impurity in a metal. 

The application of a potential V g through an electric field 
via a back gate [1] shifts the chemical potential /i and moves 
the magnetic state of the impurity in the vertical direction (y) 
in Fig. 2. We assume that x = DA/U does not change much 
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Figure 3: n T (/i) .n^fi) and x(m) for \e \/D = 0.029 and V/D = 
0.14 ( D ~ 7 eV). Left panels: x = 11 (dashed curves), and x = 5 
(solid). The impurity magnetizes inside the bubble (n-f ^ nj,). The 
vertical line marks the position of the level, eo = 0.2 eV. On the 
right: e = 0.2 eV (solid) and e = -0.2 eV (dashed) at x = 0.45. 



with applied voltage, even with screening coming from a fi- 
nite Fermi energy. Hence, the magnetization of the impurity 
can in principle be turned on and off, depending only on the 
gate voltage applied to graphene. This is better illustrated 
by looking at the behavior of the impurity magnetic suscep- 
tibility. In the presence of a field, the energy of the impurity 
spin states changes to e CT = eo — u^bB + Un- a . In the 
zero field limit, the magnetic susceptibility of the impurity, 

X = Ms J2a a ( c '"cr/di?) B= o (A 1 b is the Bohr magneton, and 
B is an applied magnetic field), can be calculated straightfor- 
wardly from Eq. dTOb : 



2 \ ^ An a 
-Vb A. 



1 - u 



de- 



l-U 2 



dn-v dn„ 



(14) 



In the lower panels of Fig. 3 we show xii 1 ) f° r e o = 0.2 
eV, V ~ 1 eV and D ~ 7 eV, for different values of x. The 
corresponding magnetization line for this set of parameters 
is defined by the solid curve with black circles in Fig. 2a. 
While the impurity remains non -magnetic for any y at x = 11 
(U ~ 40 meV), as shown in Fig. 2a, the impurity state al- 
ready crosses the magnetic boundary twice for x = 5, where 
U is nearly twice larger. In this case, a large magnetic mo- 
ment of ~ 0.5 fis forms below the energy of the level, at 
fi - 0.18 eV (Fig. 3a and 3b). At x = 0.45 (U = 1 eV), 
the local moment exists for very large // ~ 1 eV, and a strong 
and uniform magnetic moment of ~ 0.9 /is forms in almost 
the whole magnetic region (see Fig. 3c). A similar qualita- 
tive behavior for the magnetization occurs when eo < (Fig. 
3c,d). As U becomes large (> 1 eV), the magnetic transi- 
tion becomes very sharp. For eo ~ 0.5 eV and V = 1 eV, 
the impurity can magnetize for U > 0.1 eV. While the local 
charging energy U for transition metals in a metallic matrix is 
of the order of ~ 5— 10 eV [ 18], in graphene, where the effec- 
tive hybridization can be large due to the linear increase of the 
DOS, the critical U for magnetization of the impurity can be 
much smaller. Hence, transition elements and molecules that 
usually do not magnetize when introduced in ordinary metals, 
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teraction is purely ferromagnetic due to the vanishing of the 
Fermi wavevector, kp = h/vf B24I1 . However, at finite bias 
voltage the RKKY interactions display 2kp oscillations de- 
caying like 1/r 3 [25] that can couple the magnetic moments 
ferromagnetically or anti-ferromagnetically depending on the 
position and geometry of the adatom lattice (that can be con- 
veniently chosen using a STM). Hence, by changing the bias 
voltage a variety of different magnetic phases can emerge. 



Figure 4: Spectral function (in units of 1/eV) of the /-electrons at 
the Fermi energy y, for \eo\/D = 0.029 and V/D = 0.14 (D ~ 7 
eV). (a) x = 11 (solid curve) and x = 5 (dashed) for eo > (see 
Fig. 2). (b) x = 0.45 for e > (solid) and e < (dashed). 



can actually become magnetic in graphene lfl6j[l7ll . 

In order to show that the spectroscopic functions of the 
magnetic impurities can also be controlled by electric field 
effect we show, in Fig. 4, the spectral function of the local- 
ized electrons calculated at the Fermi energy: Atf(uj — /i) — 
2ir J2 a P//,o-(m)- The solid line in Fig. 4a is a non-magnetic 
resonance in a situation where the impurity state does not 
cross the magnetic boundary of the scaling diagram by chang- 
ing y (/i) for some fixed x. In the other curves of Fig.4, the 
spectral weigth splits between two peaks located around the 
magnetic transitions, near fi ~ eo and eo + U (see Fig. 3). 

The dependence of the impurity density of states with fi, 
and hence with gate bias, allows for the identification of the 
formation of local moments through ordinary transport mea- 
surements. For finite \i, the hybridization of the itinerant elec- 
trons with the localized level renormalizes the charge scatter- 
ing channels and hence the carrier conductivity, a = 2e 2 \y,\r, 
where t -1 is the impurity scattering rate. Second order per- 
turbation theory gives 11811 : 



-T 1 OC71 V 2 A//(/a), 



(15) 



where t^ 1 is the scattering rate of the electrons in the absence 
of impurities and no is the impurity concentration. In the limit 
of very large U, however, the scattering is dominated by the 
spin flip channels in the Kondo regimeJH M E3 B . 
When eo is located in the experimental range accessible by 
the application of a gate voltage ~ —0.3 to 0.3 eV U], the 
shape of the dip in the conductivity produced by the impurity 
scattering can indicate not only the position of the energy level 
but also the presence of local magnetic moments (notice that 
the non-magnetic resonance in the spectral function is quite 
symmetric). 

In the presence of a finite density of magnetic moments a 
macroscopic magnetic state can egress due to the RKKY in- 
teraction between them. At the Dirac point (u = 0) the in- 



In conclusion, we have examined the conditions under 
which a transition metal adatom on graphene can form a local 
magnetic moment. We find that due to the anomalous broad- 
ening of the adatom local electronic states, moment forma- 
tion is much easier in graphene. Furthermore, the magnetic 
properties of adatoms can be controlled by electric field ef- 
fect allowing for the possibility of using graphene in spin- 
tronics. We thank V. M. Pereira, J. M. B. Lopes dos San- 
tos, A. Polkovnikov, and S. W. Tsai for helpful discussions. 
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